Nonthermal broadening in the conductance of double quantum dot structures 
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We study the transport properties of a double quantum dot (DQD) molecule at zero and at finite 
temperature. The properties of the zero temperature conductance depends on whether the level 
attraction between the symmetric and antisymmetric states of the DQD, produced by the coupling 
to the leads, exceeds or not the interdot tunneling. For finite temperature we find a remarkable 
nonthermal broadening effect of the conductance resonance, when the energy levels of the individual 
dots are detuned. 
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I. INTRODUCTION 

Double quantum dot (DQD) systems exhibit 
a wide range of interesting and fundamental 
physical phenomena, such as Coulomb block- 
ade oscillations of conductancei^^&I&^iSA 
the formation of "double quantum dot 

mo i ccu i c ii 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. 22. 23. 24. 25 

or the Kondo effec t 26 ' 27 i 28 ' 29 . Experimen- 
tally, DQDs have been formed in semiconductor 
heterostructures^'i 2 -^^^^^^^, in single wall 
carbon nanotube o 20 ' 21 ^ 22 ' 23 and in InAs nanowire o 24 ' 25 . 
In recent years interest in DQDs has also been driven by 
the quest for a solid-state based qubit, the elementary 
building block of a quantum computer. In the presence 
of an interdot coupling t c , coherent electrons states can 
extend over the whole DQD system, resembling therefore 
the formation of chemical bonds in molecules. This 
interdot coupling t c controls the exchange interaction 
of electron spins, assumed to be localized in each of 
the dots, and hence eventually the operation of a 
corresponding solid-state qubit system^ , as well. 

Another interesting phenomenon, seen in 
experiments 16 performed in Coulomb blockade regime 
at high source-drain voltages is that the levels of one of 
the quantum dots (QDs) can act as a low-temperature 
filter for the other QD. This means that on detuning 
the energies of the levels participating in the resonant 
tunnelling process, the width of the resonance peak in 
conductance is independent of the temperature. 

Surprisingly, a complete analytic description of non- 
thermal broadening in not currently available, despite 
the fact that this effect is alluded to in the early work of 
Ref. Q] and has been demonstrated numerically in Ref. [ll| 
for coherently coupled quantum dots at large source- 
drain voltages e\V\ 3> t c . A complete analytical theory 
of coherently coupled dots is desirable, since in all the 
above theories, the shape of the resonant peak is found 
to be Lorcntzian. Experimentally non-Lorentzian line- 
shapes are found, but these are attributed to inelastic 
scattering. It is therefore of interest to ask under what 



circumstances non-Lorentzian peaks are found in the ab- 
sence of inelastic scattering. 

In this paper we consider a DQD system with co- 
herent interdot coupling as well as with coherent cou- 
plings between the dots and leads. We assume that 
the charging energies of the dots are negligible or can 
be treated as constant shifts 3 - and therefore the system 
can be described by an effective single particle model. 
(For studies of the electronic correlations in DQDs see 
Refs. laii^EiElifl) We note that experimental results 
on DQD molecules are often explained with such sim- 
ple effective single particle model o 15 ' 18 i 19 ' 20 , which usu- 
ally assume that couplings of the dots to the leads can 
be neglected. We first show how the finite interdot and 
dot-lead couplings affect the lineshape of the zero bias 
conductance resonance. Our results are non-perturbative 
and take into account all orders of interdot and dot-lead 
tunnel processes. We also consider the finite temperature 
conductance and present results for the temperature de- 
pendence of the peak height of the conductance resonance 
as well as its broadening. We show that in our model the 
non-thermal broadening effect in resonant transport can 
be observed even in the zero-bias limit, if the dot-lead 
coupling strength exceeds the interdot coupling t c . 

The paper is organised as follows. In Sec. II we intro- 
duce a very general description of coupled dots and derive 
the zero temperature transmission formula. In Sec. Ill 
we discuss the properties of the zero-temperature and 
zero-bias conductance as functions of the energy levels of 
the two dots and of the various couplings in the systems. 
In Sec. IV we derive a finite temperature conductance 
formula for the DQD system and discuss temperature- 
dependent transport properties of the system. 



II. MODEL 

We consider two quantum dots coupled to left and right 
leads with tunnel coupling Tr, and Tr, respectively as 
shown in Fig. [TJa). The interdot tunnelling coupling is 
denoted by t c and it is assumed that only one energy 



2 



level in each dot is relevant. The intradot as well as the 
interdot Coulomb interactions are neglected. 




FIG. 1: Double quantum dot coupled to left (L) and right (R) 
leads with interdot coupling t c . 



where e = (s i + i 2 )/2, T = (T L + T R )/2 is the total 
coupling strength to the leads, while Ae = (e~i — £2)/%, 
AT = (Tl — T_r)/2 are the asymmetries of the dot ener- 
gies and the couplings to the leads, respectively. 

Having obtained the transmission function T dd (E), the 
linear conductance at finite temperature is given by 



Gt 



2c 2 
h 



T dd (E) 



dfo(E) 
dE 



AE, 



(3) 



As shown in appendix [XJ the electron transmission co- 
efficient T dd (E) of the double dot system can be written 



T dd {E) = 



E-E+\ 2 \E-E_\ 2 ' 



(1) 



Here E± are the poles of the transmission: 

E± =e+iT± ^(Ae + iAT) 2 +t 2 c , (2) 



where f (E) = [1 + cxp((E - n)fk B is the equilib- 

rium Fermi distribution, /.i being the chemical potential 
of the leads. For k B T <C Ep we can take /i w Ef and ex- 
tend the lower bound of the integration to —00 in Eq. ([3]). 
The resulting integral then can be calculated by contour 
integration. The zero temperature transmission T dd (E) 
has only simple poles if |Ar| ^ t c , Ae ^ and the finite 
temperature conductance reads: 



Gt = 



2e 2 T L T R t 2 c 



h ir{k B T) A \ {uj + - cj_) 



(1) 
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where ip^(z) is the first polygamma function^, lj± = 
(E± - E F )/k B T (here E± are the poles of T dd (E) in the 
upper half complex plain, given by Eq. ©), 9 denotes 
the imaginary part, and * stands for complex conjuga- 



tion. In case of Ae = and |Ar| = t c , the transmission 
T dd (E) has second order poles and therefore Gt is given 
by 



Gt = 



e 2 T L T R t 2 c 1 
h irk B T (3w) 3 



(1) 
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Here uj = (e — £j? + iT)/k B T and V'' 2 ' ) ( 2: ) is the second 
polygamma functional. 



III. ZERO TEMPERATURE LINEAR 
CONDUCTANCE 



In what follows we briefly discuss the properties of the 
zero temperature conductance with emphasis on the ef- 
fect of couplings r^, r R and t c . The understanding of the 
zero temperature case then helps us to interpret the finite 
temperature behaviour of the conductance in Sec. IV . 



From Eq. [3J the linear conductance Go at zero temper- 
ature is given by the Landauer formulae 32 



Go 



2e 2 



T dd (Ep). 



(6) 



Using Eq. ([I]) and (|2|) one can obtain an explicit expres- 
sion for Go a s a function of Ep and the characteristic 
energies of the DQD system: E\, £2, T, Ar and t c . For 
Ar = this expression agrees with the result obtained 
in Ref. [H for serially connected dots, but otherwise also 
describes the case when Ar 7^ 0. In an experiment Ep 
would be kept fixed and the energy levels of the dots as 
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well as the tunnellings to the leads would be changed by 
side and top gates. Therefore in what follows without 
loss of generality we can set Ep = 0. 

Our aim is to understand the properties of the con- 
ductance in the £%) plane [or equivalently in (e, Ae) 
plane, where e and Ae are defined after Eq. pj]. It fol- 
lows from Eq. Jl} that depending on the ratio t 2 /A, where 
A = \J (Y 2 L + Y'j t )/2, the conductance has either one or 
two maxima in the (ei, £2) plane. As shown in Fig. [21 for 
t c 3> A the conductance is enhanced in two boomerang 
shaped regions, while for i c «A there is one maximum 
in the conductance at i\ = £2 = Ep [see Fig. [3]. 

To see the dependence of the conductance on the en- 
ergy levels of the dots and on the various couplings in the 
system, we now consider certain directions in the (ii, £2) 
plane and study the cross sections of the conductance 
along these directions. Let us first assume that the levels 
£1, £2 of the two dots are kept aligned, i.e. Ae = and 
we study Go along the e axis. 

If t c > A, upon varying e two resonances occur in the 
conductance [see Figs. [H Ufa)] at energies 



£ ± = iv/^A 1 , (7) 

corresponding to symmetric (-) and antisymmetric (+) 
molecular states. Note that in an isolated double dot 
molecule where Tl — Yr — 0, the energies of the 
symmetric (antisymmetric) states are e± = ±t c . Thus 
Eq. (0 shows that connection to the leads produces level 

attraction. Since A = \Jy 2 + AT 2 , the magnitude of 



linewidth T: 

gr-, w 2e 2 tj Y L Y R 

h (t\ - (r 2 /4 + Ar 2 )) [(£ ± y/% - Ar 2 ) 2 + r 2 ] ' 

. (8) 

Regarding the enhancement of the conductance in a 
boomerang-shape areas in the {ii^i-z) plane (see Fig. [2]), 
it is easy to prove that if Ar = 0, for a given value 
of Ae, the transmission Tdd has maximum at s = 
±-J Ae 2 + i 2 — r 2 . This is the equation of a hyperbola 
in the (e, Ae) plane and also helps to understand the ob- 
served structure of the conductance when t c 3> | Ar| 7^ 0, 
which is the case in Fig. [5] 

The separation between the two resonances decreases 
as A is increased while keeping t c fixed. Finally, the two 
resonances merge when A = t c meaning that due to the 
coupling to the leads, the energies of the symmetric and 
antisymmetric states become degenerate at this value of 
A. For A > t c the conductance then has only one res- 



10 1 , , , 1 — 1 




-10 1 1 1 J 1 — 1 

-10 -5 5 10 

Ej/f 




-10 -5 5 10 

FIG. 2: Zero temperature conductance at E = Ef as a func- 
tion of ei, £2 for t c 2> A. For details see text. 

the level attraction depends both on the total coupling 
strength Y and on the asymmetry Ar of the couplings. 
Similar level attraction has been reported in R ef. |33l for 
QDs attached in parallel to the leads and in Ref. [34] for an 
Aharonov-Bohm ring with a quantum dot in each of its 
arms. It can also be shown that Go reaches the quantum 
limit 2e 2 /h at energies given by Eq. |7]) only if Ar = 0, 
i.e. for symmetric couplings to the leads. For t c ^> A, 
when the two resonances are well separated, the lineshape 
around e — ±\Jt 2 — A 2 is approximately a Lorentzian of 



FIG. 3: Zero temperature conductance at E — Ef as a func- 
tion of £1, £2 for t c <SC A. For details see text. 

onance at energy e = Ef [see Figs. El H[b)], an d from 
Eqs. Q and has the form 

G = ^LT R t 2 c 

0[) M(e + t ) 2 + r+][(e-t~c) a + r?.]' 

Here t c = ^t 2 c - AY 2 , T± = Y if |Ar| < t c < Y and 
one can show that the conductance on resonance Go(z = 
Ep) is smaller than 2e 2 /h except for the special case 
t c = V^lTr— [see Fig. gjb)]. On the othe r hand, if 
t c < \AY\ < f, we have t c = and r± = r± ^Ar 2 - t 2 
in Eq. ©, and we find that G (e) < 2e 2 /h for all e. Thus 
for t c < A the lineshape of the resonance can again be 
approximated by a Lorentzian around e = Ep but for 
larger |e| it decreases as ~ 1/e 4 . 

Further understanding of the properties of the con- 
ductance can be gained by considering Ae 7^ 0, i.e. finite 
detuning between the levels of the (isolated) dots. In 
principle one could follow any path in the (ei,^) plane 
to study the effect of finite Ae, but let us consider two 
simple yet important cases. Let us first assume that e is 
kept fixed at the value where the conductance is largest 
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the conductance as a function of Ae reads 
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FIG. 4: (a) Conductance for t c 3> A along the e axis using 
Eq. (Tjrj (solid line). The approximation given by Eq. (fSj is 
shown with dashed lines. We used t c /Y — 7 and AT /T = 0.5. 
(b) The transmission for t c = \/T L T R < A given by Eq. 
The parameters are t c /T = 0.954 and AT/T = 0.3. For details 
see text. 



and we vary only Ae, ie we study the conductance paral- 
lel to the Ae axis. Analytical progress can be made most 
easily for Ar = 0, i.e. r L = T R = T. 

When t c > A = T and e = e± i.e. e equals the energy 
of the symmetric (antisymmetric) state, the conductance 
as a function of Ae reads 



So(Ae) 



2e 2 



AT 2 ti 



h Ae 2 (Ae 2 +4r 2 ) +4r 2 i 2 ' 



(10) 



(This corresponds to taking the cross section of T^E) 
along the line denoted by "a" in Fig. [5]) We see that the 
lineshape is basically a Lorentzian for Ae <C T, but falls 
more rapidly for Ae > T. The full width at half max- 
imum (FWHM) is Ae w 
simplifies to Ae^ w 2^/2T(t c 



2 v /2r( x /r 2 + t 2 c -r), which 



r) for t c > r. 

Another obvious way to study the effect of finite Ae 
is to fix one of the energies, (e.g. £2) and vary only s\. 
In contrast with the previous example, this means that 
both e and Ae are being varied. For Ar = 0, keeping 
e~2 = e± fixed and varying only i\ (see line " b" in Fig. \2§ 
yields a particularly simple result for the conductance: 



So(Ae) 



2e 2 



h Ae 2 +T 2 ' 



(11) 



We see that the lineshape in this case is a simple 
Lorentzian which is, interestingly, independent of t c . The 
FWHM reads Ae^ = 2I\ 

On the other hand, if t c < A = T, one can easily show 
that for fixed e = Ep the conductance along the Ae axis 
reads 



Go (Ae) 



2e 2 



AT 2 t 



2-1-2 



h (Ae 2 +T 2 +t 2 ) 2 ' 



(12) 



meaning that the lineshape is approximatelly a 
Lorentzian for Ae < ^JT 2 + t 2 and the FWHM is 



Ae w = 2^{V2- l)(r 2 + t 2 ). Finally, if we keep e 2 
aligned with Ep and vary only i\ (see line "a" in Fig. [3]) 



Go(Ae) 



2c 2 



Ae 2 



H 

4 



1 



(13) 



Thus the lineshape is a Lorentzian which means a weaker 
Ae dependence of the conductance than in Eq. fT2")) . This 
is readily seen in Fig. as high conductance ridges ex- 
tending further along the i\ and e 2 direction than along 
the Ae axis. The width of the resonance is Aew = 
{T 2 + t 2 c )/T. 

To complete our analysis of the zero temperature 
physics of our model, we now briefly discuss the local 
density of states (LDOS) in each of the dots. The LDOS 
of dot 1 (2) can be obtained from the diagonal matrix el- 
ements of the Green's function Gdd of the DQD system 
(see Eq.lA"5j): 

pi{E) = --Im{G DD ) ll: i = l,2. (14) 
Our main findings arc summarized in FigJSl If the levels 
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FIG. 5: (a) Local DOS pi(E) for t c /T = 2, Ar = (solid 
line), and pi(E), p 2 {E) for Ar/r = 0.6 (dashed and dotted 
lines, respectively), (b) Local DOS pi(E), p2{E) (solid and 
dashed lines, respectively) for AT/F = 0.6, t c /T = 0.1, and 
(c) for tc/T = 0.5_ (d) Pl (E) for Ae/T = 0.5, t c /T = 0.1 
(solid line) and t c /T = 2.0 (dashed line). 

of the two dots are equal e~i = e 2 = e then depending 
on the ratio of the coupling t c and the asymmetry of the 
dot-lead couplings Ar, one can discern three cases. For 
Ar = the LDOS, which is the same in both dots, is 
a superposition of two Lorentzians centered on energies 
E = e ± t c (see Fig. [S^a)). Upon decreasing the ratio 
t c /Ar the twin peak structure of the LDOS remains as 
long as t c / AT > 1 with the LDOS of the dot with smaller 
dot-lead coupling being is larger than that of the other 
dot's. If the coupling between the dots is weak so that 
i c /Ar < 1 [as in Fig.^b)] then the LDOS in each of the 
dots is basically a Lorentzian of width approximately that 
of the corresponding couplings T^r^). In the interme- 
diate regime of i c /Ar < 1 an interesting difference in the 
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LDOS of the two dots can be observed (Fig. E^c)). The 
dot having stronger coupling to the corresponding lead 
has a double peaked LDOS while the other's LDOS is sin- 
gle peaked. If the levels of the dots are detuned (Ae 7^ 0) 
we consider the case when Ar = 0. The LDOS then de- 
pends on the ratio t c /Ae. If the detuning of the levels 
is large compared to the coupling i.e. when t c /Ae <C 1 
the LDOS in dot 1 (2) is peaked at E = i% (e^) and has 
little overlap with the LDOS of the other dot. In the op- 
posite limit of t c / Ae ^> 1 the LDOS is a superposition of 
two Lorentzians peaked at the energies of the molecular 
states E = e ± t c . 



IV. FINITE TEMPERATURE LINEAR 
CONDUCTANCE 
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FIG. 6: Conductance (in units of 2e 2 /h) as a function of 
e for three different temperatures: ksT/T — 1,2.5,5 (blue, 
green and red lines, respectively). Curves with solid lines 
show the exact result of Eq. @ , the curves with dashed lines 
are calculated using Eq. (|17[) . 



We now consider the properties of the conductance 
for finite temperature. Although the formulas given in 
Eqs. and (|5|) are not easily readable, insight into their 
physical content can again be gained by scanning along 
certain direction in the (£1,22) plane, as in the previ- 
ous section. Let us focus on the case when |AF| =^ t c , 
Ae 7^ so that the finite temperature conductance is 
given by Eq. (HJ) and let us first address the question of 
the temperature dependence of conductance in case of 
t c < A, Ae = 0. As before, we assume that Ep is kept 
constant and e is varied. The energy and temperature 
dependence of Q F as given by Eq. Q is shown by solid 
lines in Fig. [5J As the temperature increases, the reso- 
nance gradually broadens and its height decreases mono- 
tonically with T . This behaviour can be understood in 
the limit of k B T > A > t c (with k B T < E F ). In this 
case around e — Ep, where the conductance is signifi- 
cant, one finds that \u±\/2ir = \E± - E F \/2Trk B T < 1 
and therefore the expansion of the polygamma functions 
around 1/2 can be used: 



m ,l iz s 7T 2 iipW(l/2)z tt 4 z 

W H — I ) « h— 

v K 2 2%> 2 2% 8 



■+C(z 3 ), |argz| < 7T. 

(15) 

Substituting this expansion into Eq. ([3]), we find that 
to leading order in T the peak height Gr,max decreases 
monotonically with the temperature: 



2e 2 



ll r 



i. 



h k B TT L + T R {T L T R + t 2 c ) 

if k B T > A > t c and Ae = 0. (16) 



A similar ~ 1/T decay can be found for the conduc- 
tance of a single dot 3 ^. However, the double dot re- 
sult differs from the result for a single dot by a factor 
of 2t 2 /(TLTji + t 2 .), i.e. the degeneracy of the symmet- 
ric and antisymmetric levels does not simply contribute a 
factor of 2 to the conductance as one could na'ivly expect, 
but also a factor of t 2 c /(T L T R + i 2 ). As shown in Fig. [51 



for k B T 3> A > t c the conductance is approximated by 



g T (e) S3 C/T,max COSh 



2k B T / 

if k B T > A >t c and Ae = (17) 

so that apart from the amplitude £/T,max the lineshape is 
the same as for a single dot. It follows from Eq. (fP7f that 
the FWHM is a linear function of the temperature with 
a slope of 2 acosh(v / 2) 8 '' 3 -. We see that for high temper- 
atures along the e axes, the conductance of a double dot 
is similar to that of a single dot with degenerate energy 
levels. 

In contrast to the above behaviour, if we keep e = Ep 
fixed or e"2 = Ep fixed and calculate the transmission 
as a function of Ae for different temperatures, we see 
that the resonance peak is not broadened by tempera- 
ture (Fig. [7]). Indeed, if we numerically calculate the 
FWHM as a function of T using Eq. ([4]) , we see in Fig. [8] 
(red curve) that in both cases after an initial increase of 
FWHM in the regime of k B T < T, the width of the reso- 
nance approaches a constant value as the temperature is 
further increased. For Ar = using the expansion shown 
in Eq. (JTSJ) we find that for fixed e = E F the FWHM is 
given by 



Ae w ,T « 2Vr 2 +i 2 



2\3/2 



(r 2 + t 2 ) 

2{k B Tf 



(18) 



therefore the high temperature value of the FWHM 
only weakly increases with temperature since (r 2 + 
i 2 ) 3/2 /(2/c B T) 2 < 1 [see Fig. IE]. The ratio of the zero and 
finite temperature FWHMs is Aew,T / Aew ~ 1-55 [see 
Aew after Eq. (fl~2"|) ]. On the other hand, if e.g. e2 = Ep 
is fixed and ei is changed, a similar calculation yields 

fp2 1 +2\3/2 

Ae w , T « 2^1^ - 1 ^~J> 2 , (19) 

thus Aew,T again shows a weak temperature dependence 
for k B T > r and compared to the T = case we see that 
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FIG. 7: Conductance (in units of 2e 2 /h) as a function of 
ei for £2 = Ef fixed using four different temperatures: 
fc s T/r = 0, 0.5, f .0, 2.0 (light blue, dark blue, brown and red 
lines, respectively). 



Ae WiT /Ae w » 2r/yT r +lf, which for r > i c gives 
Aew,T / Aew ~ 2. We find therefore that for both scenar- 
ios the FWHM increases as a function of temperature if 
fcgT < r but for k B T 3> T it approaches a constant value 
of 2-\/r 2 + if. This is a remarkable non-thermal broad- 
ening effect and is a central result of our paper: while 
the peak height decreases monotonically with tempera- 
ture [see Eq. (fT6|) ]. the width of the peak, when changing 
only Ae or §i, approaches a constant value, as shown 
in Eqs. ([IS]). ([15]) and in Fig. U We emphasize that al- 
though for the analytic calculation we assumed Ar = 0, 
our numerical results show that the nonthermal broaden- 
ing is also present for a finite difference in the couplings to 
the leads, i.e. when Ar ^ 0. One can also show that for 
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FIG. 8: FWHM as a function of temperature, along the 
Ae axis. The red curve shows the numerically obtained 
FWHM using Eq. (j4|, the green curve is the function given 
by Eq. (IS). We used t c /T = 0.5. 

large T, where the FWHM is approximately constant, the 
lineshape of the resonance is approximately a Lorentzian: 



2e 2 7iT 



r. 



h k B T Ae 2 + T 2 + t 2 



if k B T > T >t Cl and Ar = 0. (20) 
This result holds for both of the scenarios discussed i.e. 



either e or £2 being fixed while Ae is varied. 

Let us now consider the case of t c > A. The tempera- 
ture dependence of Qt as given by Eq.0]along the e axis is 
shown in Fig.[9| We see that the conductance peaks grad- 
ually broaden and their height decreases with increasing 
temperature and finally they merge into a single peak for 
k B T > t c . The maximum of the conductance can then 
be found at e = Ep and increasing the temperature fur- 
ther to the regime of k B T 3> t c > A this peak behaves 
the same way as the case where we assumed t c < A. 
This occurs because in the high temperature limit k B T 




FIG. 9: Conductance as a function of e for three different 
temperatures: ksT/F — 0, 1, 8. (red, green and blue lines, re- 
spectively). The ratio of t c and Y is t c /T = 10, while Ar = 0. 



is the largest energy scale in the system and |cj±| <^k B T. 
Hence the expansion shown in Eq. (|15[) is applicable and 
leads to the same results as Eqs. (16]), ([17]), (T8]) and (fT9]l . 

For t c 3> k B T,A however, when one can still observe 
two distinct peaks in the conductance (see e.g. the green 
curve in Fig. [9]) different approach has to be employed. 
Let us focus on the peak at the energy of the symmet- 
ric state, i.e. e_ = — \Jt 2 c — A 2 . (Analogous consider- 
ations can be made around the energy e+ = \/t 2 c — A 2 
of the antisymmetric state.) If t c ^S> k B T,A, then for 
energies around e_, one finds that |w + |/2-7r <C 1, but 
|w_|/27r ^ 1 in Eq. (U). This implies that in the expres- 
sion of Qt the polygamma functions, whose argument is 
can be neglected compared with the other two terms, 
which are functions of uj + . Indeed, from the expansion 
of polygamma functions for large arguments 



7 + 7^ + 773 
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00, |argz| < 7r, (21) 



it is clear that the contribution of •^^ 1 ^(l/2— iu>-/2ir) and 
its complex conjugate's in Eq. (0]) is very small if e s» e_ 
which means that for a zero bias measurement one can 
consider the symmetric state as a single resonant level. 
The finite temperature conductance around e_ can be 
then approximated by 



7 



Qt 

if t c > k B T, A and Ae = 0. 



2e 2 



r f ,r 



h n{k B T) 2 T L + T R ^tl - Ar 2 V (w+ - w_) 



(i) 



2tt / 



+ c.c. 



(22) 



In case fc^T is much larger than A, i.e. for 
t c S> k B T 3> A one can use the expansion shown in 
Eq. (fT"5]) to obtain further approximations of Eq. (|22|) . 
For the peak height we find very similar results to the 
case of t c < A. To leading order the peak height de- 
creases monotonically with temperature: 

t/T.max- h 2k B TT L +T R {T L T R + t 2 c ) 

if i c > k B T > A and Ae = 0. (23) 

Note that there is a factor of 1/2 difference compared to 
Eq. (fl6j) . because the symmetric state behaves as a single 
resonant level. One also finds that the lineshape is rather 
well approximated by 

if t c > fc B T > A and Ae = 0. (24) 

especially for e -C £-■ However, due to the other res- 
onant level at energy e+, the lineshape is in fact not 
symmetric around e_ as Eq. (|24|) suggests. Comparision 
of the exact result given by Eq. (01 with the approxima- 
tions of Eq. and Eq. (JH]| is shown in Fig. [TUl A 



Ae^K.T for each case can be calculated using Eq. (f2"2"]l and 
the expansion shown in Eq. (|15p. The resulting formulas 
are however rather complicated and not too informative. 
Numerical calculations shown in Fig. [11] clearly indicate 
that the resonance is broadened by temperature in the 
regime of t c 3> k B T 7 T. 



a) 0.6i 1 1 b) 35 
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FIG. 11: The conductance as a function of Ae, for three dif- 
ferent temperatures: ksT/T — 1,5,10 (blue, green and red 
lines, respectively) and fixed ej = e_ (a). The numerically 
calculated FWHM using Eq. Q as a function of temperature 
(b). We used Ar = and t c /T = 20. 

We can conclude therefore that in the studied sys- 
tem the non-thermal broadening of the conductance res- 
onance occurs only if k B T is the largest energy scale. 
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FIG. 10: Conductance as a function of e around the energy of 
the symmetric state e_ if t c 3> feT 3> T. The exact result of 
Eq. Q is shown with red line, the approximation of Eq. (I22p 
with blue and the formula given by Eq. (|24[) with green line. 
Parameters: fc s T/r = 8, t c /k B T = 5, AF L = 0. 

small but noticable deviation of the approximation given 
by Eq. (j2"4"|) from the exact result Eq. (01 can indeed be 
observed for e ^ e_, while Eq. ([22]) gives a better ap- 
proximation over the whole energy range around e_. 

For the properties of the conductance along the Ae 
axes or along the <?i axes, the characteristic energies 



V. CONCLUSIONS 

In conclusion, we have studied the linear conductance 
of a double quantum dot molecule at zero and finite tem- 
peratures. We have found that the coupling of the dots to 
the leads produces level attraction, which depends both 
an the total coupling strength Y and on the asymme- 
try Ar of the couplings. We have discussed the prop- 
erties of the conductance when the interdot coupling t c 
is larger (smaller) than this level attraction. In partic- 
ular, at zero temperature we have given explicit expres- 
sion for the line shape of the conductance in the 
£2) plane along certain experimentally important axes. 
Considering the finite temperature conductance we have 
discussed the temperature dependence and the lineshape 
of the conductance along these axes. We have showed 
that if the temperature is the largest energy scale in the 
system, the conductance resonance, which arises due to 
the detuning of the energy levels of the quantum dots, 
is not broadened by the temperature. Our results can 
be relevant for understanding of those recent experimen- 
tal results where an effective single particle description is 
adequate. 
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APPENDIX A 

In this appendix we derive the electron transmission 
coefficient through serially coupled, coherent quantum 
dots, connected to multi-channel leads. There are nu- 
merous equivalent approaches to computing transport 
through such phase-coherent structures, including recur- 
sion methods and transfer matrix technique o 42 ' 43 . Here 
we employ the Green's-function method and notation 
presented in Ref. 0, in which the Hilbert space is di- 
vided into a sub-space A containing the external leads 
and a sub-space B containing the two dots. 

We start by considering isolated left and right dots, 
which are each described by a single quantum state |/i > 
and I/2 >, with energy levels e\ and £2 respectively. 
When these are coupled together by a Hamiltonian hi2, 
the 2x2 Greens function gs of the coupled dots is given 
by 



9u 



~t* 922 



(Al) 



where t c =< /1I/112I/2 > and gj^ = E — Ej. (This rep- 
resentation is convenient, because the self-energy matrix 
in Eq. (|A5|) below is then diagonal.) 

The effect of coupling the left (right) dot to the left 
(right) lead via a coupling matrix W\ (W2) is represented 
by self energies = ol — iTl (£7? = <?r — ^Fr), (where 
<tl, Tl, gr, Tr are real) defined by 

S L = < M W !\ n L > 9A(n L ) < n L \W x \h > (A2) 



which yields the transmission coefficient Tdd via the for- 
mula 



T dd = 4Tr\r{L)G DD r(R)Gfi DD ] = 4r L r«|(G^) 12 | 2 , 
where 



T(L) 



and 



T(R) 



r L 





r R 



(A6) 
(A7) 



(A8) 



Finally from Eq. 



{G DD )i2 = W[(£-£i-a L +ir L )(£-£2-^+Tfl)-|t c | 2 ], 

(A9) 

and writing i\ — £i+ol and £2 = £ 2 +ctr, yields equation 
([1]) of the main text. 

We note that this equation resembles Eq. (20) of 
Ref. H]]. However the latter omits the self energy terms 
ul and an, which in general are non-negligible. 



APPENDIX B 

Introducing the dimensionless variable 
y = (E — E F )/k B T, the integral I in Eq. (J3J) reads 



4r L r i? t c 2 



dvf'iv) 



(k B T) 4 J_ 00 (y - uj+)(y - cj*)(y - LU-)(y - luI) 

(Bl) 

where f'(y) = l/cosh 2 (|) and uj ± = (E± - E F )/k B T. 
This integral can be calculated using contour integration. 
Care has to be taken, however, because the integrand is 
not bounded on the imaginary axes. Nevertheless, one 
can calculate this integral as a sum of two contour inte- 
grals, as shown in Fig. [T2l 



and 



St? = ]T < f2\W±\n R > g A (n R ) < n R \W 2 \f2 >, (A3) 



where \n L rm > is a channel state belonging to the L (R) 
lead and 3a("-l (r)) is the channel Green's function, such 
that 

9a (R) = \ n L(R)> 9 A { n L (R)) < n L (R)\ (A4) 

riL («h) 

is the corresponding surface Green's function. 

In the presence of the leads, the Green's function Gbb 
of the double dot is given by 



U DD — 9d 



s. 



(A5) 




FIG. 12: The two integration contours to calculate the inte- 
gral in Eq. © • Filled circles denote the second order poles of 
— 9f g^ , open circles show the poles of Tdd(E). 



Closing the contours in the upper half plane, the con- 
tributions of the two contours along the imaginary axes 
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cancel, except around the (second order) poles of the 
derivative of the Fermi function fo(E) (shown by filled 
circles in Fig. [121). These poles are located at (2k + l)z7r, 
where k is an integer. The other contribution to the inte- 
gral comes from the poles E± of the transmission function 
Tdd(E) (denoted by open circles in Fig. [T2|) . Summing 
all the contributions from the poles and using the series 
expansions of the first polygamma function 31 



V> (1) (z) = 



E 

k=0 



1 



(z + ky 



(B2) 



and of the 1/ cos 2 (z) function 



1 oo 

cos z (z) z — ' 



fe=0 



1 



1 



((2k + 1)tt - z) 2 ((2k + 1)tt + z) 2 
one can finally obtain the result shown in Eq. 
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